Statistics of Coulomb Blockade Peak Spacings within the Hartree-Fock Approximation 
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We study the effect of electronic interactions on the addition spectra and on the energy level 
distributions of two-dimensional quantum dots with weak disorder using the self-consistent Hartree- 
Fock approximation for spinless electrons. We show that the distribution of the conductance peak 
spacings is Gaussian with large fluctuations that exceed, in agreement with experiments, the mean 
level spacing of the non-interacting system. We analyze this distribution on the basis of Koop- 
mans' theorem. We show furthermore that the occupied and unoccupied Hartree-Fock levels exhibit 
Wigner-Dyson statistics. 

PACS numbers: 72.10.Fk, 73.20.Dx 



Introduction — It has recently been established that 
studying the low temperature fluctuations of the con- 
ductance through quasi-isolated nanostructures is an ex- 
cellent tool for probing electronic interactions p). In 
such small devices called artificial atoms or quantum 
dots it is possible to control and successively vary the 
total number of electrons from zero S to a few hun- 
dreds [||-0. When coupled weakly to leads, the conduc- 
tance of the dots is vanishingly small between pronounced 
peaks because of Coulomb blockade S. By tuning a 
gate voltage connected to the dot a conductance peak 
is observed whenever the ground state energy E G {n) of 
the dot containing n electrons becomes degenerate with 
the energy of the dot with n + 1 electrons. At almost 
zero drain to source voltage this situation is expressed as 
Hn = Ec(n + 1) - E G (n) = /iLoad, where ^ Lca d and n„ 
is the chemical potential of the leads and of the dot con- 
taining n electrons, respectively. The spacing between 
adjacent conductance peaks is given by the difference 



Mn- 



-P»-1 = E G (n + l) - 2E G {n) + E G (n-l) . (1) 



Recently, a new generation of experiments on Coulomb 
blockade in ballistic and diffusive quantum dots has 
shown that the conductance peaks fluctuate both, with 
respect to their heights M and spacings ^-0]. Most of 
the features of the peak height statistics are well under- 
stood within random matrix theory (RMT) |1C(| . How- 
ever, all measured |0-M] peak spacing distributions resem- 
ble a Gaussian form, while the constant interaction model 
H , which properly accounts for the average peak spacing, 
together with RMT predicts a Wigner-Dyson statistics. 
Moreover, the observed widths of the distributions vary 
between the experiments from a width comparable to the 
mean single-particle spacing A B to a width consider- 
ably larger plkl, whereas the widths of the GOE and 
GUE distributions are about A/2. These experiments 
suggest that the widths scale rather with the charging 
energy E c than with A. 



Different theoretical studies addressed this issue: RPA 
diagrammatic perturbation theory yields corrections to 
the constant interaction predictions of order 1/g pl| or 
1/sfg P4 I, where g is the dimensionless conductance. 
However, these approaches, which are valid at small r s 
(high densities), cannot easily explain the large widths 
(~ 8A) of the latest experiments @. A refined RMT 
approach [Of, which accounts for shape deformations of 
the dot while adding electrons, explains the Gaussian 
profiles, while the widths remain comparable with A. 

Numerical exact diagonalization studies of an interact- 
ing tight-binding model of ~ 10 electrons show that for 
r s = 1 the distribution is indeed Gaussian 01. This is 
also the situation when spin is included Jljjl^]. Never- 
theless, questions were raised regarding the relevance of 
these calculations to the experimental system due to the 
small number of electrons considered, the use of a tight- 
binding model for which it is difficult to relate the param- 
eters to the experimental systems, and due to the small 
values of g (of order one) accessible in this approach. 

In this paper we study interacting spinless electrons 
on continuous two-dimensional disordered cylinders at 
zero temperature. The interactions are considered within 
the self-consistent Hartree-Fock (SCHF) approximation, 
which enables us to treat up to n = 50 electrons. By 
using a continuous model, the experimental parameters 
can directly be related to the calculations, and we obtain 
higher values of g {g ~ 10) close to the experimental val- 
ues. We show that within the SCHF approximation the 
peak spacing statistics is Gaussian in the presence or ab- 
sence of magnetic field. The variance of the fluctuations 
is much larger than predicted in the independent-particle 
picture. This is in agreement with the exact diagonalisa- 
tion results, with a density functional approach [ pT[ , and 
with other recent HF studies of disordered tight-binding 
models |M and clean quantum dots [O. Using Koop- 
mans' theorem [081 we show that the calculated widths 
and Gaussian profiles can, at least partly, be related to 



specific Hartree interaction integrals (for not too strong 
interactions). We further show and explain that, on the 
contrary, the spacings among the occupied or the unoc- 
cupied HF levels follow GOE and GUE level statistics. 

Numerical procedure — We begin with the HF equa- 
tions for spinless electrons in real space, 



H + V dis (f) + -%- I 

47T6 J n , 



ELilVy(rO| 



A|2 



-dfi' 



ipk(r) 



47T6 



Y.Ur 3 {r')Mr')Ur) ^ . _ ... 

— J - J — — — d\V = e k i/j k (r) . (2) 



Here e is the dielectric constant, Ho = — {h 2 /2m* e )d 2 / dr 1 

and Vdis(j) = Eni^/f^ - Pq) the potential from N s 
(5-scatterers at random locations r q homogeneously dis- 
tributed on the cylinder surface $7. The strengths X q are 
chosen from a box distribution. 

We use < 9 < 2tt as polar and < y < n as vertical co- 
ordinates and boundary conditions, ipj(0, 0) —ijjj(9, it) — 
0, Vj(0 + 27r,y) = e i27r, ^ o ^(0,2/), to account for a flux 
4> piercing the cylinder. <fta = he/ e is the flux quantum. 
The set of Eqs. (0) is solved self-consistently by diagonal- 
izing the Hamiltonian matrix in the basis of Hq. To this 
end we expand %(f) = e~ i6<t>/(t ' il)j{f) = J2i a\ J) Ui{r) in 
the periodic basis functions Ui(f) — {l/ir)e me6 sin(n y y) 
of Ho, where ng = 0,±1,±2, ... and n y — 1,2,3... are 
good quantum numbers which define a level I — l(ng, n y ) 
of Hq- We emphasize that the enumeration of the lev- 
els l(ng,n y ) does depend on flux. Using energy units 
h /{m* e R 2 ), where R is the cylinder radius, the matrix 
elements of the Hamiltonian (J2J) read 

H(h,h) ^ elS lul2 +V dis (h,h) + uY,Ai [ ,i'Vi [ ,i' 2 (li,h) ■ 



Here, 1,1^1^ denote eigenstates of Ho and e° = 
1/2 [(rag + (j>/(f)o) 2 + (n y /2) 2 ~\ is the scaled energy of an 
electron on a clean cylinder of aspect ratio unity. 
Vdis(h,h) are the disorder matrix elements. U = 
R/2rB determines the interaction strength in dimen- 
sionless units with tb = 47re?i 2 /(m*e 2 ) the Bohr ra- 
dius in the sample material. Ay ^ — Yll=i a v a V 
and V ViA {l x ,h) = {hJ'MhA) 1 - ih,l' x W2M) with 



1/lr 



[sii 



e')/2) + (y-y') 



/\21 



-1/2 . 



in real 



space. While the matrix A changes along the iterations 
until self-consistency is reached, Vd; s is calculated only 
once. We used 450 x 450 matrices and ensured that the 
results are practically unaffected by this truncation. 

Generally, interaction effects are governed by the 
charge density n s and by U. We used U = 0.2tt 2 = 1.97. 
On the one hand, this choice still allows convergence 
of the SCHF, on the other hand this is in the experi- 
mental parameter regime prtl: In GaAs re = 10~ 2 ^m 
(to* = 0.067m e , e = 13eo). This gives for a typical dot 



area of (2nR) 2 ~ 0.07 - 0.3/^m 2 the estimate U ~ 2 - 4 
which is similar to the numerical value used. Accord- 
ingly, typical experimental values of r s ~ 1 for GaAs 
||-£| and r s ~ 2 for Si are also close to our numerical 
value r s — 2U \Jir/n that varies between 1.0 and 1.6 for 
20 < n < 50 electrons. The disorder strength used corre- 
sponds to g ~ 10 describing the experimental situation. 
g was calculated from the inverse participation ratio p9[ . 
Results and analysis — For a system containing n elec- 
trons the SCHF ground state energy is given by |20| 



n _. n 

£ G (n) = £ £fe -- Y, Uk\v\j,k)-(j,k\v\k,j)] 



(3) 



fc=i 



j,k=i 



where j and k enumerate SCHF levels and we redefined 
v = U/\f — r'\. Eq. (ra) was used to calculate the peak 
spacings s n in Eq. (n^T The peak spacing distributions 
were generated from an ensemble of 600 quantum dots 
(with different disorder realizations) and sequences of s n 
with 20 < n < 50. The results are depicted as the his- 
tograms 1 in Fig. 1 for = and Fig. 2 for finite <fi. 
The curves are centered around the mean peak spacing 
(s n ) ~ 8.3A. Clearly each distribution fits a Gaussian, 
curves 3, defined by the average and variance of the data. 
In the following we analyze the numerical results. As- 
suming that the single-particle states are unchanged by 
adding (subtracting) an electron to (from) the many- 
electron system we can approximate the difference in en- 
ergy of the ground states of the system with n + 1 and n 
electrons according to Koopmans |18J by 



lh, 



E G (n + 1) - E G (n) ~ e n+1 (n) 



(4) 



Here, e n+ i(n) is the n+lst energy eigenvalue in the sys- 
tem of n>l electrons. In view of Eq. (tih, we then have 
approximately 



£ n +i(n) -e„(n- 1) 



(5) 



The numerical distributions of the s n from Eq. (H) are 
depicted as curves 2 in Figs. 1 and 2. They show con- 
siderable agreement with the corresponding distributions 
(curves 1) of the SCHF ground state energies, Eq. p), 
and their Gaussian fits. This shows that, for the interac- 
tion strengths considered, a peak spacing analysis based 
on Koopmans' relation is justified. As visualized in the 
right inset of Fig. 1, the peak spacing is related to en- 
ergy differences from systems with different number of 
electrons pl| . We emphasize that e„(n — 1) ^= s n (n). 
Their difference is typically large (50% and more) on the 
scale of s n in Eq. (H). Hence s n ^ e n +i(n) — s n (n). 

We also checked directly the validity of Koopmans' re- 
lation (0) by comparing the two numerical values for \i n 
and 6,1+1 (n). The typical relative error, the difference 
(e n +i(n) — fj, n )/ ' Hn, is always positive and ~ 0.4%, indi- 
cating the validity of Eq. (0). 




FIG. 1. Distributions of normalized Coulomb blockade 
peak spacings x = s n /A for weakly disordered cylindrical 
quantum dots and flux <j) = 0. The histograms 1 and 2 are cal- 
culated from the self-consistent Hartree-Fock ground state en- 
ergies (Eq. Inh) and by using Koopmans' relation dq). Curve 3 
is a Gaussian defined by the average (the charging energy) and 
variance of the data of curve 1. Curve 4 is the distribution of 
the exceptional direct integral (g) yielding a considerable con- 
tribution to s„ in Eq. ([]). Curve 5 is a Gaussian fit. Left in- 
set: The nearest neighbour spacing distribution of SCHF sin- 
gle-particle energies for the occupied and unoccupied states is 
compared with the GOE statistics p(s) — (tt/2)s exp(— 7ts 2 /4) 
(dotted line). Right inset: A typical realization of SCHF 
eigenvalues. The dotted line connects eigenvalues relevant for 
calculating e.g., s n =4,2 using Koopmans' relation (J5J). 




FIG. 2. SCHF peak spacing distributions as in Fig. 1 but 
for 4>/<j>o — 0.25 which reduces the fluctuations. Inset: The 
nearest neighbour distribution, as in the left inset of Fig. 1, 
of SCHF single-particle energies is compared with the GUE 
distribution p(s) = (32/tt 2 )s 2 exp(— 4s 2 /-7r) (dotted line). 

To understand the origin of the fluctuations of s„ we 
recall that the HF energies Sk depend on n and obey |2C| l 

e k (n) = (k\H +V di8 \k)+J2U,k\v\j,k)-(j,k\v\k,j). (6) 



Using this relation in Eq. (||) we obtain 

s n ~ [(n + 1\H + V dis \n + 1) - {n\H + V dis \n)] 
+ (n, n + l|i;|n, n + 1) — (n, n + l\v\n + 1, n) 

n-l 

+ ^2 [U, n + l\v\j, n + 1) - (j, h\v\j, n)] 

3=1 

n-l 

- Yl [U^ + l\v\n + l,j)- (j>n\v\n,j)] . (7) 
i=i 

The tilde sign denotes states calculated in the system 
with n — 1 electrons. When v — > one recovers the 
Wigner-Dyson distribution for s n which is modified upon 
increasing the interaction. The two sums in Eq. (R) con- 
tain differences between direct terms and between ex- 
change terms which partly cancel each other. This is not 
the case for the direct and exchange integral in the second 
line which therefore represents an exceptional contribu- 
tion to s n . The direct term dominates over the exchange 
term because the integrand of the former is always pos- 
itive. Hence we expect that the exceptional direct in- 
tegral, (n,n + l\v\n,n + 1), describing the interaction 
between the n-th and n-l- 1st electron and given by 



U 



|*„(r)| 2 |vI/„ +1 (r-)p 



■dridn' , 



(8) 



o Jn 1 



significantly contributes to s n . 

It has been established 22 23] that for moderate disor- 



der the eigenfunctions fluctuate in space with small cor- 
relations. They even survive for different eigenfunctions 
but are weaker p3j. This should hold also in our SCHF 
calculations because the screening of the disorder poten- 
tial is not perfect and disorder effects will persist. Be- 
cause of the long-ranged Coulomb term, the central limit 
theorem should hold fairly good for the direct (Eq. (0)) 
and exchange integrals (despite the small wave functions 
correlations) and one thus expects to find a Gaussian 
distribution. For the interaction strengths used (U ~ 2) 
the direct term (g) can dominate the s n which implies a 
Gaussian distribution for the peak spacings, too. 

In Figs. 1 and 2 the numerical distribution of the di- 
rect term m) is shown as curve 4 and is compared with 
a Gaussian fit, curve 5. This distribution is expected 
to be centered to the left of the full numerical spacing 
distribution, because we have also to account for a shift 
on the scale of the mean HF single-particle spacing A^f 
according to the first line in Eq. (R). Note that below 
the highest occupied and above the lowest unoccupied 
HF level, the spacing is A HF = 0.314 > A = 0.183. The 
shift to the left in Fig. 1 and Fig. 2 is - 0.280 - 0.88A ffF 
and ~ 0.270 ~ 0.84A#f, respectively. 

Due to the appreciable interaction strength the discus- 
sion based on Eq.(J7T) is practically independent of flux. 



Both s n statistics are close to Gaussians but with flux- 
dependent widths. 

Comparison with experiment — The Gaussian form of 
the distributions is in good agreement with all experi- 
mental results [Q-0. The average peak spacing of the 
interacting systems is (s„)/A = 1.52/0.183 = 8.27 > 1 
(curve 1 in Figs. 1 and 2), similar to the experimental 
situation in GaAs |J-|(|. (s n ) is also in agreement with 
the estimated (scaled) classical capacitance of the cylin- 
der, [e 2 /4ne(TrR)]/[h 2 /m* e R 2 ] = 0.4tt = 1.26 

The width (RMS) of curve 1 in Fig. 1 (Fig. 2) equals 
1.6A = 0.19(s„) (1.3A = 0.16(s„)). This is considerably 
larger than the widths 0.0966 = 0.52A (0.0743 = 0.41A) 
for the non-interacting case which are in good agreement 
with GOE (GUE) predictions ||. 

The enhancement of the widths is compatible with the 
first GaAs experiments pLpfl and in line with the Si ex- 
periment M] which shows a larger width for higher values 
of r s . The width is large compared to the experiment 
H where it was close to its non-interacting value. Since 
all experimental settings are quite similar, it is not clear 
what can lead to the different widths observed in IJ . 

SCHF single-particle statistics — We finally consider 
the statistics of the HF level spacings Sk with k ^ n to em- 
phasize that they are not related to the Gaussian spacing 
distribution. By using Eq. (o) we have 

s k = e k +i(n) -£fe(n) 

= (k + 1|F + V dis \k + 1) - (k\H + V dis \k) 

n 

+ ^2[(j,k + l\v\j,k + l)- (j,k\v\j,k)] 

3=1 

n 

-^2[(j,k + l\v\k + l,j)-(j,k\v\k,j)]. (9) 

3=1 

In contrast to Eq. (m) the exceptional and dominant 
direct integral does not appear in Eq. (m leading to 
St < s„. It is easy to see that for k > n all terms con- 
tribute and for k < n all terms except those with j = k 
and j = k + 1 that cancel each other. Hence there is no 
essential difference in the s k distributions in the two k re- 
gions. Because of the tendency of the Coulomb terms to 
cancel each other, the terms in the second line in Eq. (H) 
are expected to dominate the average. The left insets of 
Figs. 1 and 2 demonstrate that the numerical distribu- 
tions of the normalized s k follow indeed GOE and GUE 
statistics (dotted lines), respectively, showing the effect 
of time-reversal symmetry breaking. This RMT type of 
statistics is in line with other HF results that was found 
in a different type of experiment [|26| . 

Conclusion — We have studied Coulomb blockade 
peak spacing fluctuations for interacting spinlcss elec- 
trons within the self-consistent Hartree-Fock approxima- 
tion. The main features of the conductance peak spacing 
fluctuations are similar to the experimental ones. The 



Gaussian distribution of the addition spectrum is inter- 
preted as the result of the dominance of direct interaction 
terms governed by spatial fluctuations of the eigenfunc- 
tions. On the other hand the Hartree-Fock single-particle 
levels were shown to follow Wigner-Dyson statistics. 
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